Astronomy & Astrophysics manuscript no. a2218cosmo 


February 2, 2008 


(DOI: will be inserted by hand later) 





Multiple-Images in the Cluster Lens Abell 2218: 
Constraining the Geometry of the Universe ? 
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Abstract. In this Letter we present a detailed study of the lensing configuration in the cluster Abell 2218. Four multiple- 
images systems with measured spectroscopic redshifts have been identified in this cluster. These multiple images are very 
useful to constrain accurately the mass distribution in the cluster core, but they are also sensitive to the value of the geometrical 
cosmological parameters of the Universe. Using a simplified maximum likelihood analysis we find < C1 M < 0.30 assuming 
a flat Universe, and < Am < 0.33 and w < - 0.85 for a flat Universe with dark energy. Interestingly, an Einstein-de Sitter 
model is excluded at more than Act. These constraints are consistent with the current constraints derived with CMB anisotropies 
or supernovae studies. The proposed method constitutes an independent test of the geometrical cosmological parameters of the 
Universe and we discuss the limits of this method and this particular application to Abell 2218. Application of this method 
with more sophisticated tools and to a larger number of clusters or with more multiple images constraints, will put stringent 
constraints on the geometrical cosmological parameters. 
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1. Introduction 

The present Cosmology framework is characterized by a num- 
ber of parameters which sets the global geometry of the 
Universe, its history and dynamics. The quest for these pa- 
rameters is a long-standing issue in Observational Cosmology 
and is still the main driver of a large number of experiments. 
Combining constraints coming from the power spectrum of the 
CMB anisotropies and the luminosity distances of distant type 
la superno vae (SNIa), a new s tandard model of cosmology is 
emerging llSpergel et alJl2003b : a flat Universe with an accel- 
erating expansion (Qm - 0.27 and Qa 0.73). To quantita- 
tively explain these results the concept of dark energy has been 
put forward, characterized by the ratio of pressure and energy 
density w = Px/pxc 2 , which reduces to the vacuum energy 
(the cosmological constant) for w — - 1 . T here is however, no 
strong observational constraints on w yet JSpereel et alj F2003 
give only w < - 0.78). 

This new standard cosmology is getting very popular. 
Although the flatness of the Universe see ms robust, the ex- 
act value of Qm is still a matter of debate (iBridle et alJ l2003: 
iBlanchard etalTl2003h as it is essentially driv en by the SNIa 
results which can be discussed jRowan-Robinsonll2002l) . In 
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order to independently probe the large scale geometry of the 
Universe, we propose to explore the potential use of cluster 
lenses as a long range optical test bench. Preliminary analy- 
sis of this method was first detailed by iLink & Piercd {1 998) 
using simple lens models. Recently, we extended their work 
using more detaile d simulations of realistic clusters of galaxies 
dCiolse et alJl2002h . The basic idea of this method is that each 
set of multiple-images identified in a cluster lens strongly con- 
strains the cluster potential. As the scaling of the mass model 
depends on the ratio of the angular distances £>ls/A}s, it will 
also depends on the geometrical cosmological parameters (Qm, 
Qa and w). In order to constrain these parameters, the com- 
bination of several sets of multiple images in a single lens is 
mandatory to disentangle between the degeneracies in the lens 
model. With a minimum of 4 systems of multip le images with 
known spectroscopic redshifts, we showed JGolse et alj E002) 
that one can put reasonable constraints in the (Qm, Qa) plane 
with some charact eristic degeneracies in the fitted parameters 
JGolse et all2002h . 

In this Letter we apply this lensing test to the well stud- 
ied cluster-lens Abell 2218. Section 2 describes details of the 
lens modeling and the cluster mass distribution. The results of 
the optimization are discussed in Section 3 and a conclusion is 
presented in Section 4. When necessary we scale the physical 
parameters with Q M = 0.3, Qa = 0.7, h = 0.65 with Hubble 
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Table 1. Redshift and properties of the multiple images sys- 
tems in Abell 2218. Nj is the number of images detected in 
the HST image and used in the optimization (note that a third 
image is predicted for the system #4 but is too faint to be de- 
tected). Na = 2(Nt - 1) is the number of constraint s intro - 
du ced in the optimizat ion. References: (1) | Pell6 et al 1 ill 992b . 
l^ lEbbels et all i ll 9961) . G^ lEllis et alJJ200lh . 



Multiple image system 


zs 


N t 




References 


#1 


0.702 


4 
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(1) 


#2 


1.034 


3 


4 


(1) 


#3 


2.515 


3 


4 


(2) 


#4 


5.576 


2 


2 


(3) 



constant Ho = 100 h km s 1 Mpc 1 . Thus at the cluster redshift 
Z = 0.176, 1" corresponds to a linear scale of 2.08 h^ 1 kpc. 

2. Lens Modeling 

2.1. Leasing and other constraints 

Abel l 22 18 is one of the richest clusters in the Abell catalog 
lAbell etalJll98 9|). A s pectroscopic survey of the galaxies by 
iLe Borgne et alJ Jl992) led to an average redshift zl = 0.1756 
and a galaxy velocity dispersion <x = 1370*™ km s _1 (based 
on 50 cluster members). It is also one of the few clusters for 
which both an accurate lens model and the identification of 4 
systems of lensed multiple images with spectroscopic redshift 
are presently available. Figure 1 and Table 1 show this list of 
multiple images with their properties. A number of strong lens- 
ing m odels have been discussed in the lite rature (Kneib et al. 
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ll996HAllenll998tlEllis et al 12001 UNataraian et all2002hh us- 
ing only part of the current available constraints. 

Comparisons of the inferred lensing mass with the mass 
distribution derived from the X-ray emission of the cluster 
pointed to a strong discrepancy between th e two estimators 
(Miralda-Escude & Babull ll995l lAllenlll998B . if one assumes 
hydrostatic equilibrium. A recent and detailed study of the gas 
distri bution in this cluster us ing high quality Chandra observa- 
tions (Machac ek et al.l l2002) partly confirms this discrepancy, 
especially near the center. It demonstrates that A 22 18 is prob- 
ably not fully relaxed. However a much better agreement is 
found in the outer parts of the clust er where the X-ray mass 
compares with weak lensing masses ( Squi res et al.ll996tlAllenl 
1998). Our analysis is however independent of the physical 
state of the intra-cluster gas. 

2.2. The different gravitational lens components 

We start from the model described in Knei b et all {l996): 
the cluster mass is distributed within two halos of dark mat- 
ter centered respectively on the main cD galaxy and on 
the second brightest galaxy. The mass profile of each halo 
is modeled with a so-called truncated PIEM D ( "Pseudo- 
Isothermal Elliptical Mass Distribution", Kassiola & Kovner 



I1993L iKneib etalj fl996) characterized by 7 parameters: 4 are 
geometrical (center (Xq, Yq), ellipticity e and position angle 6) 
and 3 describe the mass profile (velocity dispersion cr , core 
radius r c and truncation radius r t ). Note that our definition of 
the ellipticity is s = (a 2 - b 2 )/(a 2 + b 2 ). The numerical sim- 
ulations of many different lens configurations and their fits by 
several analytical mass distributions have shown that cosmo- 
logical constraints are not very sensitive to the choice of the 
analytic model: provided there is enough freedom in the num- 
ber of parameters, the fits can easily adjust the true mass profile 
at the location of the multiple images. For illustration if a char- 
acteristic core radius is included among the free parameters, it 
will shrink to very small values if the t rue profile is sing ular 
(see some examples in the simulations of lGolse et alJ2002l) . 

Furthermore, we include the contribution of the 37 bright- 
est cluster galaxies with magnitude m < 19.5 (i.e. m»+2) and 
we associated a truncated PIEMD to each of them. The geomet- 
rical parameters (center, ellipticity, position angle) are fixed to 
those of galaxy light parameters, and the mass profile parame- 
ters are scal ed with the total mag nitude, using the prescription 
proposed by Nataraian & Knei ri ll 19971) and inspired from the 
standard Faber-Jackson and Kormendy relations: 
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<To*, G c * and 8 t * are reference values for am* = 17.5 galaxy. 
For the last relation, a = 0.5 means that a constant M/L ratio 
applies to all galaxies, while for a + 0.5, M/L scales as L a ~ l t 2 . 
To minimize the number of parameters in the model, we fix C * 
to a very small value 8„ = 0.048" (or a physical scale r„ = 
0. l/i 1 kpc) for a nearly singular mass profile. Indeed, galaxies 
are essentially characterized by their central velocity dispe rsion 
and the extension of their halo (Nataraian & Kneib 1997). 

2.3. Optimization of the main lens parameters 

In order to avoid biases in the determination of the cosmologi- 
cal parameters, we use the following procedure, al ready te sted 
and e valuated by previous numerical simulations JGolse et alJ 
2002): using a sparse sampling in the (Qm, A) plane (with steps 
of 0. 1 for both parameters), we optimize all the other model pa- 
rameters with a Monte-Carlo method. Thus we do not bias the 
final optimization toward a given cosmology. For each of the 
2 main dark matter halos, we adjust the 7 parameters of the 
PIEMD and for the individual galaxies we explore the 3 pa- 
rameters: o"o„, #t* an d a - 

A wide range of values is allowed for each parameter dur- 
ing the Monte-Carlo initialization. When scanning the area 
[0 < Om < 1,0 < Q.a < FJ, the minimum^ 2 is found for 
Cl M = 0.001 and Q A = 0.9, with a reduced y 2 . = 6.05. 

m /v ' A. mm 

Similarly, assuming a flat Universe and scanning the area [0 < 
Qm < 1; -1 < w < 0] the minimum is located at Q. M - 0.101 
and w = -1, with a reduced y 2 =6.65. These high y 2 . 

' A mm to t min 

values are likely representative of yet non-perfect mass models 
and possibly also of underestimation of intrinsic errors, espe- 
cially in the image positions. Anyhow, in both cases, the values 
of the parameters which describe the two halo potentials and 
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the galaxies are close to those obtained hv lKneih et al l Jl99(t> 
and lNataraian et all J2002bl) . The main halo is centered on the 
central cD galaxy with a shift of a few arc seconds with respect 
to it. The second halo is well positioned on the second bright- 
est galaxy, again with a few arc seconds shift with respect to its 
center (Fig.^. 

The orientations of the two main halos are slightly differ- 
ent from those of the light distribution but they show a clear 
tendency of alignment between two halos. Moreover, the ellip- 
ticity of the main halo (e\ 0.28) is significantly smaller than 
that of the cD isophotes (e c o = 0.563) while the ellipticity of 
the secondary halo is quite l arge (e? . 61). T hese features 
were already pointed out by iKneib et all Jl996h in their lens 
model. They may represent the signature of a merging phase 
of the secondary halo the main mass concentration, in good 
agreement with the recent X-ray analysis of the Chandra data 
bv lMachacek et alJ J2002I ). 

The characteristic values of galaxy halo mass distribu- 
tion are compatible with current values deduce d from galaxy- 
galaxy lensing analysis like those found by iNataraian et alJ 
(2002a) in their study of 6 clusters of galaxies. The a pa- 
rameter which represents the variation of the truncature radius 



with luminosity is about 0.9. The value 0.5 which corresponds 
to a constant MI L value for all galaxies is ruled out, as al- 
ready claimed by Nataraian et alJ J2002bh and the M/L ratio 
then scales as L 0A . The brightest galaxies seem to have a more 
extended and massive halo than fainter ones, as already sug- 
gested in studies of the fundamental plane o f elliptical galaxies 
^Dressier et alii 987U.Torgensen et al Jl 99ft. 



2.4. Cosmological constraints 

Starting from the lens model determined in the previous sec- 
tion, we perform a second level of optimization within the 
(£2m,£2a) or me (Qm,w) plane, with steps of 0.01. Only the 
most critical parameters are optimized in this modeling, the 
others being fixed at their previously determined "best value". 
The fitted parameters are the velocity dispersions cr 01 and <Xo2, 
and the core radii 9 C [ and 6 C 2 of the two main halos. When 
changing the cosmology, we keep the lens efficiency of the 
other clumps identical by rescaling the velocity dispersion of 
the galaxies so that ctq oc (cr i + 0-02), therefore reducing the 
number of free parameters. 
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Fig. 2. Confidence levels obtained in the (Qm, Qa) plane (top) 
and in the (Qm, w) plane (bottom) for a chi-square computed 
with v = 10 degrees of freedom (see text for details). Contour 
levels vary from lcr to 4<x from the darkest to the lightest lev- 
els, i.e. probability levels are 68.3%, 95.5%, 99.7%,and 99.99% 
respectively. 



The results of this optimization are displayed in Fig. [5] 
Only the best model and its x 2 value are kept at each step, 
with no real marginalization on the lens parameters. This is 
an approximation which is justified by analogy with the anal- 
ysis of CMB data: in the case of Gau ssian distributions of the 
errors, Tegm ark & Zaldarriaeal (|2000) demonstrate the equiv- 
alence between a full multidimensional marginalization and a 
much simpler maximization. In our case, it is not clear whether 
the probability distribution function of the lens parameters is 
close to Gaussian, but to save computer time, we decided to 
use a similar approach. Further analysis of the full likelihood 
distribution is necessary to fully validate the method. But is is 
out of the scope of this paper, presented mainly as a demon- 
stration case applied on a single cluster lens. In Fig. 13 the con- 
fidence levels are given for a number of degrees of freedom 



Table 2. Estimates of the uncertainty in the lens parameters de- 
rived from the lcr x 2 contours in the cosmological parameters 
space. 





cr i (kms ') 


a-Q2 (km s ') 


0ol (") 


0c2 (") 


(Qm, £2a) 


1039!^ 
1035+}^ 


383*29 
382+f 


18.60!^ 
18.60-- 


8.05!™° 
7 qo+i.m 

7 -0.06 



v = Nc - N L = 16 - (4 + 2) = 10, where Nc is the number 
of constraints displayed in Tableland Nl is the number of pa- 
rameters in the optimization (4 for the cluster model and the 
2 cosmological parameters). The contours of these confidence 
levels follo w very clear l y the expected degeneracy except at 
low Qa (see lGolse et alJ 120021) for a detailed discussion about 
these degeneracies). Therefore we are confident a posteriori 
that the modeling of the cluster is a fair representation of the 
true potential and mass distribution. 

In addition, if we overlay on the lcr x 1 contours the con- 
tours of the fitted lens parameters (croi, o"oa, Oci, &ci), we derive 
an estimate of the variation of these parameters for the "best 
models". Table 2 summarizes these values considered as an es- 
timate of the error bars of the lens parameters. 

The consequences on the constraints of the cosmological 
parameters are encouraging. It is clear from this analysis that an 
Einstein de Sitter Universe is excluded, at more than 5 cr. And 
as expected, the constraints are more stringent on Qm than on 
Qa- If we assume a flat Universe we get some narrow windows 
on the parameters: 

Q M < 0.22 or Q A > 0.78 

in close agreement with the constraints derived from 
Supernovae experiments. And if we assume the existence of 
a dark energy component in a flat Universe, we find 

Q M < 0.37 and w < -0.80 

These results are comparable to those obtained by c ombining 
CMB and supernovae data (w < -0.78. ISpergel et alJ J2003h ). 
Our results also compare well with other recent determinations 
issued from the statistics of gravitational lenses. For example, 
the CLASS (Cosmic Lens All Sky Survey) survey of radio 
galaxies provided constraints on Qm and Qa or w ve ry simi- 
lar to ours, but with larger error bars JChae et alJ2002l) . 

However some limitations in the procedure are quite ob- 
vious. Sources of uncertainties are shared between uncertain- 
ties in the image positions (even with the accuracy of Hubble 
Space Telescope (HST) images), and errors in the mass models. 
Indeed, although the dependence in the mass profile has been 
addressed in our previous paper, with accompanying simula- 
tions, the reality of the mass distribution in clusters of galaxies 
is likely to be more complex ISand et alJ2003l) . Further exami- 
nation of these limitations are in progress and will be presented 
in a forthcoming dedicated paper. 
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3. Conclusion and future prospects 

We have shown in this paper that we can derive reasonable 
cosmological constraints from the very detailed analysis of the 
lensing configuration of the cluster of galaxies A2218. The nec- 
essary conditions for this study are simple: deep multicolor 
HST images of a well selected cluster-lens, identification of 
a minimum of 4 families of multiple-images systems and se- 
cure redshift measurement of each family, which ranges from 
Z =0.702 to 5.576. With these constraints, and provided the 
mass distribution can be modeled by the sum of a dominant 
component and smaller additional ones (all following a trun- 
cated PIEMD mass profile), the geometrical problem can then 
be solved. The cosmological constraints presented in this pa- 
per are of similar accuracy than those derived from Supernovae 
analysis. Interestingly, both analysis are purely geometrical and 
completely independent tests, but they are not sensitive to the 
same combination of distances thus providing nearly orthogo- 
nal constraints in the (Qm, ^a) plane. 

This new method to constrain the cosmological parameters 
is very attractive, especially in view of the outstanding perfor- 
mances of the Advanced Camera for Surveys (ACS) on board 
of HST and the development of Integral Field spectrography 
that allow to secure the redshift of many multiple images in a 
ve ry efficient way. The very spectacular ACS images presented 
bv lBenitez et alJ (E002) on Abell 1689 show that in a very near 
future, we can use the proposed method as a very serious cos- 
mological test, by focusing on those clusters with more than 4 
multiple images with spectroscopic redshift. One advantage of 
this method is its relatively low-cost in terms of telescope time 
and relatively easy to implement - although progress is needed 
to thoroughly explore the parameter space of the mass models 
and to implement a fully comprehensive likelihood analysis. 
Such improvements are currently under investigation and will 
in a near future allow a better treatment of this exciting prob- 
lem. 
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